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A Geant4 hadronic process verification suite has been designed to test and optimize Geant4 hadronic models 
in the cascade energy range. It focuses on quantities relevant to the LHC radiation environment and spallation 
source targets. The general structure of the suite is presented, including the user interface, stages of verification, 
management of experimental data, event generation, and comparison of results to data. Verification results for 
the newly released Binary cascade and Bertini cascade models are presented. 



1. INTRODUCTION 

The Geant4 toolkit [| includes a collection of mod- 
els and packages for hadronic physics which are ap- 
plicable to various particle transport problems Q. In 
Geant4, final state generation is separated from the 
access and use of cross sections and from tracking. It is 
therefore possible and desirable to have independent, 
alternative physics models. In order to optimize these 
models and to identify their most appropriate applica- 
tion, detailed comparisons between experimental data 
and model predictions are required. 

The cascade energy region from the reaction thresh- 
old to incident hadron energies of a few GeV is prob- 
lematic for all existing Monte Carlo packages. In order 
to verify and further develop Geant4 models in this en- 
ergy range, a verification suite has been created. The 
stages of the verification, including experimental data 
handling, event generation and comparison to data 
are described here. Results of the comparison are also 
presented. 



2. VERIFICATION SUITE 
2.1. Method 

The verification suite is generic, as it is based on 
an abstract interface to a final state hadronic interac- 
tion generator. This interface can be found in Level 
2 of the Geant4 hadronic physics framework It 
focuses on quantities relevant to the LHC radiation 
environment and spallation source targets. 

The modular structure of Geant4 allows the gener- 
ation of single events with a known incident particle 
energy and an explicitly defined hadronic final state 
generator. The kinematics of secondaries produced in 
the interaction are analyzed and the resulting angu- 
lar, momentum, energy, and baryon number spectra 
are stored in histograms. The energy-momentum bal- 
ance can be controlled as well. The histograms are 



compared to published measurements of the differ- 
ential and double differential cross sections, da/dE, 
da/dtt, d 2 (7 1 'dEdfl, and the invariant cross sections, 
Ed 3 <r/d 3 p. 

The cross section contributing to the i-th bin of the 
histogram are given by 

ACT, = (TtotNi/N, (1) 

where (Jtot is the total cross section for the interaction 
being tested, N is the number of simulated events 
in the sample, and iVj is the number of times bin i 
is incremented during the run. Each bin represents 
a small region of phase space such as Ail, AE, or 
AQ.AE, into which the secondary particle goes after 
being produced in the interaction. The double diffcr- 
cntion cross section is cstimatatcd as 

dEdn AEAQ' 1 ' 

The verification suite is organized into a number of 
test cases, each defined by a unique incident beam en- 
ergy and a single target nucleus. All data files, macro 
files, kumac files and results for each case are stored 
in a separate subdirectory. The AIDA abstract inter- 
face |4| is used in the suite for histogram creation and 
filling. The output is done in ASCII and HBOOK 
formats. 

2.2. Cross Sections 

The verification is done by comparing simulation 
results with experimental data mainly from the EX- 
FOR database 5] . Only data with absolute measure- 
ments of the differential cross sections are utilized in 
the suite. The data downloaded from the database are 
re-formatted in order to provide intermediate files ac- 
ceptable for PAW analysis. In some cases re-binning 
of the data was performed during this process. This 
was required mainly for the low energy part of the 
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Table I Neutron production by incident protons. 



Target nucleus 


Beam energy(MeV) 


Be 


113, 256, 585, 800 


C 


113, 256, 590 


Al 


22, 39, 90, 113, 160, 256, 585, 800 


Fe 


22, 65, 113, 256, 597, 800 


Ni 


585 


Zr 


22, 35, 50, 90, 120, 160, 256, 800 


Pb 


35, 65, 113, 120, 160, 256, 597, 800 



Table II Pion production by incident protons. 



Target nucleus 


Beam energy(MeV) 


H 


585 


D 


585 


Be 


585 


C 


590 


Al 


585, 730, 1000 


Cu 


730 


Ni 


200, 585 


Pb 


585, 730 



spectra. If, within a given test case, only double dif- 
ferential cross section data are available, the single 
differential cross sections are obtained by numerical 
integration. 

Initial verifications have been performed for neu- 
trons and pions produced by protons incident upon 
various targets. For these test cases the initial proton 
energies were all below 1 GeV. In this energy region 
the inclusive reaction channel 

p + A -»• n + X, (3) 

has been studied experimentally for many years. The 
secondary neutrons can be identified and their ener- 
gies can be measured with good precision by using 
time-of-light techniques. Because this reaction is im- 
portant for many applications, a significant number of 
test cases have been created for it (Table 0. 

For HEP applications and in particular for the LHC 
detector, simulating secondary pion production is im- 
portant. Model verification for pion production by 
incident protons is available for several test cases (Ta- 
ble [nj for the reactions 

p + A -» n ± + X. (4) 



2.3. User Interface 

The user interface is implemented by macro files 
which allow the specification of various parameters of 
the verification. These include: 



• initial particle and its energy, 

• beam energy spread, 

• target nucleus, 

• hadronic interaction generator, 

• generator options, and 

• histogram types and bins. 

Both linear and logarithmic binning are available. 
Each verification test case contains a script which exe- 
cutes the macro files. Hence the user can perform the 
complete verification process by issuing a single com- 
mand. The results are stored in HBOOK format and 
can be processed in PAW by prepared kumac files. 

3. HADRONIC MODELS 

The verification suite has been used extensively dur- 
ing the release phase of the new Geant4 Bertini cas- 
cade and Binary cascade packages Q- The Bertini 
cascade model is a classical cascade code which was 
described in detail at another presentation at this con- 
ference 0. 

The Binary cascade introduces a new approach to 
cascade calculations. The interaction is modeled ex- 
clusively on binary scattering between reaction par- 
ticipants and nucleons. The nucleus is described by 
a detailed 3-dimcnsional model in which the nucleons 
are explicitly positioned in phase space. Free hadron- 
hadron elastic and reaction cross sections are used to 
define the collisions. Propagation of the particles in 
the nuclear field is done by numerically solving the the 
equation of motion. 

The cascade begins with a projectile and a descrip- 
tion of the nucleus, and terminates when the both 
average and maximum energy of all particles within 
the nuclear boundary are below a given threshold. 
The remaining nuclear fragment is treated by pre- 
equilibrium decay and de-excitation models 0. 

4. RESULTS OF VERIFICATION 
4.1. Inclusive neutron spectra 

The single and the double differential inclusive neu- 
tron spectra are very sensitive to the physics model 
used in the cascade code. As shown in Fig. ^ the 
Binary cascade reproduces the data rather well for all 
targets, except for energies below 50 MeV where neu- 
tron evaporation is important. In contrast, the Bertini 
cascade (Fig[2J does well below 50 MeV for all but the 
lightest nuclei. Above 50 MeV only the general trend 
of the data is reproduced, while the overall normal- 
ization improves as A increases. 
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Figure 1: Neutron yield produced by 256 MeV protons. 
Histograms - Binary Cascade predictions, points - data 
i 



Figure 2: Neutron yield produced by 256 MeV protons. 
Histograms - Bertini Cascade predictions, points - data 



Fig. [5]shows forward scattering in which the proton 
transfers its energy to one target neutron. Here the 
Binary cascade describes the data above 50 MeV. The 
Bertini cascade (FigsEJ does less well, significantly 
underestimating the data at higher neutron energies. 
This discrepancy is reduced for heavier nuclei. 

Additional verification test cases indicate that the 
Binary cascade approach is reasonably accurate for 
other angles and energies as well. Fig. [5] shows com- 
parisons for aluminum at forward and backward neu- 
tron angles at energies of 113, 256, 585 and 800 MeV. 
At 113 and 256 MeV, agreement at all angles is good. 
At 585 and 800 MeV the backward angle spectra are 
not well-reproduced. The same set of plots is shown in 
Fig. [5]for iron, and in Fig. 0for lead. For both targets 
the same trends apply as were observed in aluminum. 
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4.2. Pion production 

Another useful test of the cascade codes is to look 
at the double differential cross sections for the produc- 
tion of 7r + and 7r~. Experimental cross sections show 
that 7r + production by protons is significantly larger 
than that for ir~ . This feature is well-reproduced by 
the Binary cascade as shown in Fig|S] The trend of 
the cross section versus energy is also reproduced, al- 
though the overall normalization is underestimated by 
a factor of 2-3 for carbon, aluminum and nickel. 

Similar plots for tt + and n~ are shown for an inci- 
dent proton energy of 730 MeV. FigslUl and 1101 show 
7r~ and 7r + at forward and backward angles from alu- 
minum, while Figs. llll and [T2l show 7r~ and n + from 
copper. In all these cases the observed ratio of tt + 
to 7r~ production is reproduced. As was the case for 
597 MeV incident protons (Fig. |5J, the trend of cross 



Figure 3: Double differential cross-section for neutrons 
produced at 7.5 degrees by 256 MeV protons. Histograms 
- Binary Cascade predictions, points - data 0. 

section versus energy is also reproduced. However at 
730 MeV it is seen that the Binary cascade does an 
increasingly poor job of reproducing the overall nor- 
malization as the pion angle increases. 

4.3. Future verification work 

While the existing verification suite has already 
been used extensively, further development is re- 
quired. So far the validity of the cascade codes has 
only been demonstrated for incident protons at the 
low end of their energy ranges. In order to test the 
full range of energies and particles, many more test 
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Figure 4: Double differential cross-section for neutrons 
produced at 7.5 degrees by 256 MeV protons. Histograms 
- Bertini Cascade predictions, points - data 



cases are required, including: 

• incident proton energies up to 15 GeV, 

• incident pion energies up to 15 GeV. 

• incident neutrons, and 

• proton inclusive spectra. 



5. CONCLUSIONS 

A verification suite for hadronic interaction models 
in the cascade energy region has been designed and 
implemented. This suite has been used both to de- 
velop the Bertini and Binary cascade model codes for 
the Geant4 toolkit and to compare their predictions 
with data from thin target scattering experiments. 

The suite has proved to be an efficient tool for 
model validation which can be used for the devel- 
opment and testing of other hadronic interaction 
models. The suite is also being expanded to include 
tests for more types of incident particles and higher 
energies. 
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Figure 5: Double differential cross-section for neutrons produced in proton scattering off aluminum. Histograms 
Binary Cascade predictions, points - data S Fit I Hill if . 
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Figure 6: Double differential cross-section for neutrons produced in proton scattering off iron. Histograms - Binary 
Cascade predictions, points - data [si l9lllcllTl| . 
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Figure 7: Double differential cross-section for neutrons produced in proton scattering off lead. Histograms - Binary 
Cascade predictions, points - data fsl l9lHcllTl| . 
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Figure 8: Double differential cross-section for pions produced at 45° in 597 MeV proton scattering off various materials. 
Histograms - Binary Cascade predictions, points - data [T^|. 
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Figure 9: Double differential cross-section for ir produced in 730 MeV proton scattering off aluminum. Histograms - 
Binary Cascade predictions, points - data |T^I . 
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Figure 10: Double differential cross-section for n produced in 730 MeV proton scattering off aluminum. Histograms 
Binary Cascade predictions, points - data |T^|. 
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Figure 11: Double differential cross-section for n produced in 730 MeV proton scattering off copper. Histograms - 
Binary Cascade predictions, points - data |T^I . 





Figure 12: Double differential cross-section for 7r produced in 730 MeV proton scattering off copper. Histograms - 
Binary Cascade predictions, points - data |T^|. 
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Figure 13: Double differential cross-section for n produced in 730 MeV proton scattering off aluminum. Histograms 
Binary Cascade predictions, points - data |T^I . 
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Figure 14: Double differential cross-section for tv + produced in 730 MeV proton scattering off aluminum. Histograms 
Binary Cascade predictions, points - data |T^I . 
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